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Abstract 

We present ab initio calculations which show that the direct-band-gap, effective masses and 
Fermi velocities of charge carriers in ZnO monolayer (ML-ZnO) in graphene-like honeycomb 
structure are all tunable by application of in-plane homogeneous biaxial strain. Within our 
simulated strain limit of ±10%, the band gap remains direct and shows a strong non- linear 
variation with strain. Moreover, the average Fermi velocity of electrons in unstrained ML-ZnO is 
of the same order of magnitude as that in graphene. The results promise potential applications of 
ML-ZnO in mechatronics/straintronics and other nano-devices such as the nano-electromechanical 
systems (NEMS) and nano-optomechanical systems (NOMS). 
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1 Introduction 

Nanomaterials whose electronic properties 
can be controlled by mechanical strain are 
highly desirable for applications in nano- 
electromechanical systems (NEMS). If such 
materials happen to have direct band gaps, 
then those can find applications in nano- 
optomechanical systems (NOMS). If mechanical 
strain strongly affects the electronic properties 
of a material, then that material can be used as 
a strain sensor. Bulk ZnO, in its most stable 
wurtzite structure at ambient pressure, has 
many useful properties [TJ [2] such as direct 
and wide band gap (Eg = 3.37 eV at room 
temperature and 3.44 eV at low temperatures), 
large exciton binding energy (60 meV), strong 
piezoelectric and pyroelectric properties, strong 
luminescence in the green-white region, large 
non-linear optical behavior, high thermal con- 
ductivity, radiation hardness, biocompatibility, 
and so on. These properties make ZnO an 
excellent candidate for a variety of applications 
in optics, electronics and photonics, sensors, 
transducers and actuators. Several ZnO nanos- 
tructures in the form of very thin nanosheets, 
nanowires, nanotubes and nanobelts have been 
synthesized and characterized using different 
preparation methods [3, 4J. In 2006, using the 
density functional theory (DFT) calculations, 
Freeman et al. [5] first predicted that when 
the layer number of ZnO films with (0001) 
orientation is small, the wurtzite structures are 
less stable than a phase based on 2D ZnO sheets 
with a layer ordering akin to that of hexagonal 
BN. In 2007, Tusche et al. [B] reported the ob- 
servation of 2 monolayer (ML) thick ZnO (0001) 
films grown on Ag(lll) by using surface X-ray 
diffraction and scanning tunneling microscopy. 
Very recently graphene-like honeycomb struc- 



tures of ZnO have been prepared [7] on Pd(lll) 
substrate. Theoretical studies [5.8HT0] have now 
established the stability of the monolayer and 
few-layers (FL) of ZnO in planar honeycomb 
structures. This stability is attributed to the 
strong in-plane sp 2 hybridized bonds between 
the Zn and O atoms. 

Recent theoretical studies on two dimensional 
(2D) flat ZnO nanostructures reported (i) the 
elastic, piezoelectric, electronic, and optical 
properties of ZnO monolayer (ML-ZnO) |10j . 
(ii) room-temperature half-metallic ferromag- 
netism in the half-fluorinated ML-ZnO |11| . 
(hi) fluorination induced half metallicity in 
few ZnO layers [12], (iv) fully-fluorinated and 
semi-fluorinated ZnO sheets |13j and (v) strain- 
induced semiconducting-metallic transition for 
ZnO zigzag nanoribbons [2]. Here, we report 
our DFT based investigation of the effect of 
homogeneous in-plane biaxial strain on the 
electronic properties of ML-ZnO - a specific 
study not found in the existing literature. 
This study simulates an experimental situation 
where ML-ZnO is supported on an ideal flat 
stretchable/flexible substrate. 

2 Computational Methods 

The calculations have been performed by us- 
ing the DFT based full-potential (linearized) 
augmented plane wave plus local orbital (FP- 
(L)APW+lo) method [15J which is a descendant 
of FP-LAPW method [16]. We use the elk- 
code [17] and the Perdew-Zunger variant of LDA 
|18j . the accuracy of which have been success- 
fully tested in our previous works on graphene 
and silicene |19j (silicon analog of graphene), ger- 
manene [20] (germanium analog of graphene) and 
graphene/h-BN heterobilayer [21]. The plane 



1 



wave cut-off of |G + k| ma:r = 9.0/R m t (a.u. -1 ) 
{R m t = the smallest muffin-tin radius) was cho- 
sen for plane wave expansion in the interstitial 
region. The Monkhorst-Pack [22] k-point grid 
size of 30 x 30 x 1 was used for all calcula- 
tions. The total energy was converged within 
2/ieV/atom. We simulate the 2D-hexagonal 
structure of ML-ZnO as a 3D-hexagonal super- 
cell with a large value of c-parameter (= |c| =40 
a.u.). The application of in-plane homogeneous 
biaxial strain 5 up to ±10% was simulated by 
varying the in-plane lattice parameter a(= |a| = 
|b|);<5 = (a — ao)/ao ; where ao is the ground 
state in-plane lattice constant. Figure 1 depicts 
the top-down view of a ML-ZnO in planar con- 
figuration. 




Figure 1: Top-down view of a ZnO monolayer in 
graphene-like honeycomb structure in ball-stick 
model. The large (pink) ball represents Zn atom 
and the small (green) ball represents the O atom. 



3 Results and Discussions 

Our calculated LDA value of ao = 3.20 A cor- 
responds to the Zn-0 bond length dzn-o = 
ao/\/3 = 1.848 A, which is in agreement with 
the reported experimental [6] value of dz n -o 
1 .92 A and theoretical values of dz n ~o = 1-86 A 
by Tu and Hu [8], 1.895 A by Topsakal et al. [9], 
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Figure 2: Bands and total DOS of unstrained 
ML-ZnO within LDA. Ei? is the Fermi energy. 

1.853 A by Tu [ID], 1-85 A by Wang [12J. Our 
3.75% underestimation of dz n -o value with re- 
spect to the experimental value is due to the well 
known problem of underestimation of the lattice 
constant within LDA. In our previous studies 
|18pi9j of graphene, silicene and germanene, we 
have demonstrated that the value of c-parameter 
chosen in the construction of supercell for simu- 
lation of 2D hexagonal structures also affects the 
value of in-plane lattice constant: larger value of 
c-parameter yields a slight smaller value of a. 
Since we use a different method and a different 
value of c-parameter, the disagreement of our re- 
sult on ao with other theoretical result (!) is ac- 
ceptable. 

The electronic band structure and total den- 
sity of states (DOS) plots of ML-ZnO are de- 
picted in Figure 2. As seen in Figure 2, ML- 
ZnO is a direct band gap (E s = 1.68 eV, LDA 
value) semiconductor with both valence band 
maximum (VBM) and conduction band mini- 
mum located at the T point of the hexagonal 
Brillouin Zone (BZ). However, the actual band 
gap is expected to be larger as LDA is known 
to underestimate the gap. Our calculated LDA 
band gap of 1.68 eV is in agreement with pre- 
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Figure 3: Variation of the direct band gap of ML- 
ZnO with in-plane homogeneous biaxial strain 5. 



vious calculations of 1.68 eV by Topsakal et 
al. [9], 1.762 eV by Tu [TO], 1.84 eV by Wang 
et al. [13J. Using the GW approximations, re- 
cently Tu [TO] estimated the direct band gap of 
ML-ZnO at 3.576 eV, which showed ML-ZnO 
as a wide band gap semiconductor. Although 
both LDA and generalized gradient approxima- 
tions (GGA) do not yield the band gap correctly 
these are powerful enough to predict the correct 
trend in variation of the gap |23H26| . Since we 
focus on the trend as well as the nature of vari- 
ation in band gap (such as the possible tran- 
sition from direct to indirect gap-phase as re- 
ported theoretically for uniaxially strained ZnO 
nanotubes [27] and ZnO nanowires and nan- 
otubes [28]) rather than its absolute value, we 
employed the computationally simpler and less 
time-consuming LDA for investigating the effect 
of biaxial strain on band parameters of ML-ZnO. 

Our calculated results on homogeneous biax- 
ial strain-induced modifications of the band gap 
of ML-ZnO are depicted in Figure 3, which shows 
a strong non-linear variation of band gap with bi- 
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Figure 4: Variation of effective masses of the 
electron m* (el) , light hole m* (lh) and heavy hole 
m*(hh) at the V point of the BZ with in-plane 
homogeneous biaxial strain. 
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Figure 5: Variation of Fermi velocities of the 
electron -up (el), light hole t»F(lh) and heavy hole 
VF(hh) at the T point of the BZ with in-plane 
homogeneous biaxial strain. 
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axial strain. Within our simulated strain limits 
(« -10% < 5 < 10%) as shown in Figure 3, the 
band gap E 9 remains direct and its value varies 
with S (= [(a - a )/a )] x 100 %) as 

E g (S) = B + B 1 x5 + B 2 x5 2 + B 3 x5 3 (1) 

where B = 1.68195 eV, Si = 0.00105 eV, 
B 2 = -0.00336 eV, B 3 = 8.57346 x 10~ 5 eV, ob- 
tained by polynomial fit of our computed data. 
This may be useful in calculating the LDA band 
gap corresponding to a particular value of biax- 
ial strain 5 in the given range. Our prediction of 
strain-engineered band gap of ML-ZnO, if ver- 
ified by experiments, may be useful in fabrica- 
tion of NEMS and NOMS based on a graphene- 
analogue of ZnO. The energy bands very close 
to the T— point are parabolic, which enabled us 
to define and estimate the effective masses. Fig- 
ure 4 depicts the nature of variation of effective 
masses of electrons m*(el), light holes m*(lh) and 
heavy holes m*(hh) at the T-point of the BZ 
(see Figure 2) with our simulated biaxial stain. 
Slightly away from the parabolic bands close to 
the T— point, the bands are linear in k which 
allowed us to calculate the Fermi velocity (vp) 
values from their slopes. The nature of varia- 
tion of our calculated Fermi velocities of elec- 
trons v F (el), light holes ^(lh) and the heavy 
holes v F (hh) near the T-point of the BZ with 
our simulated biaxial stain is depicted in Figure 
5. We found that the Fermi velocity v F (el) of 
unstrained ML-ZnO (v F (el) = 0.836 x 10 6 m/s 
for 5 = in Figure 5) near the T-point of BZ 
is close to he experimentally observed value of 
v F = 0.79 x 10 6 m/s in monolayer graphene 
(MLG) deposited on graphite substrate [29] and 
the theoretically calculated DFT value of v F = 
0.833 x 10 6 m/s in MLG [30] and v F ~ 0.8 x 10 6 
m/s in MLG supported on monolayer of hexago- 
nal boron nitride [21,31]. Since the charge carri- 



ers enter into ballistic transport limit when they 
have high velocities in the absence of scatter- 
ing [32], we expect ballistic transport in ML-ZnO 
due to the high velocity of charge carriers anal- 
ogous to the case of graphene. 

4 Conclusions 

Our all-electron full-potential density functional 
theory based calculations show that band- gap, 
effective masses and Fermi velocities of charge 
carriers in ML-ZnO in graphene-like honeycomb 
structure are mechanically tunable by applica- 
tion of in-plane homogeneous biaxial strain. Our 
predictions on strong non-linear variation of di- 
rect band gap with biaxial strain (Figure 3), 
non- linear increase of m*(el), and almost lin- 
ear decrease in m*(lh) and m*(hh) with increas- 
ing strain (Figure 4) and almost linear decrease 
in the Fermi velocities of charge carriers with 
increasing strain (Figure 5) are not only novel 
but also experimentally testable by application 
of biaxial strain to a flat stretchable/flexible 
substrate supporting ML-ZnO. Our results are 
significant in that an analogous honeycomb- 
structure 'graphene' does not show any vari- 
ations of its band gap under similar biaxial 
strains and currently a lot research is being 
done not only to open but also to control the 
band gap in graphene. Owing to its direct band 
gap with strong non-linear variation with biax- 
ial strain, ML-ZnO should have potential appli- 
cations in mechatronics/straintronics and other 
nano-devices such as the NEMS and NOMS. 
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